Time-dependent resonant tunneling for a parallel-coupled double quantum dots 
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We derive the quantum rate equations for an Aharonov-Bohm interferometer with two vertically 
coupled quantum dots embedded in each of two arms by means of the nonequilibrium Green's func- 
tion in the sequential tunneling regime. Basing on these equations, we investigate time-dependent 
resonant tunneling under a small amplitude irradiation and find that the resonant photon-assisted 
tunneling peaks in photocurrent demonstrate a combination behavior of Fano and Lorentzian reso- 
nances due to the interference effect between the two pathways in this parallel configuration, which 
is controllable by threading the magnetic flux inside this device. 
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I. INTRODUCTION 

The investigation of quantum coherence in mesoscopic 
systems has been the subject of considerable interest in 
the solid state physics during the last years. In the in- 
terference experiments with a quantum dot (QD) em- 
bedded in one arm of the Aharonov-Bohm (AB) ring, 
the tunneling through a QD was proved to be coherent 
by detecting the flux-periodic current oscillations i More 
recently, HoUeitner et al. extended this idea to measure 
the AB oscillations of the mesoscopic ring containing two 
coupled QDs inserted in each of the two armsi^ Further- 
more, this parallel-coupled QDs structure has been in- 
vestigated in the Kondo regime, and an observation of 
the transition between different quantum states has been 
reported.-^ii. Clearly, the possibility to manipulate each 
of the QDs separately and the application of the mag- 
netic flux provide more controllable parameters for de- 
signing the transport properties. This has been discussed 
by several theoretical works for the stationary transport 
by nonequilibrium Green's function (GF)i^4SiLSiS 

On the other hand, time-dependent tunneling through 
coupled QDs in series has received large attention both 
theoretically and experimentally. A theoretical study of 
the photon-assisted tunneling (PAT) in double QDs given 
by Stoof and Nazaro\iiS and Hazelzet et alm^, based on 
the quantum rate equation approach, predicted that the 
photoresponse of the system exhibits satellite resonance 
peaks due to PAT processes which involved the emission 
or absorption of one photon to match the energy differ- 
ence between the discrete states of the two QDs. W. G. 
van der Wicl et al}^' measured the PAT current through 
weakly coupled QDs and discovered clearly the predicted 
extra resonance peaks under irradiation of microwave. 
Motivated by this perfect match between the theory and 
experiments, we intended, in this paper, to study the 
PAT in the parallel-coupled QDs by the quantum rate 
equations approach. In this conflguration, the additional 
bridges between the QDs and leads allow the electron 
wavefunction propagating along different pathways, then 
lead to the interference effect between them, which is dis- 



played by the fundamental AB oscillation in the presence 
of a magnetic field. Therefore, the central point of our 
study is to explore how the interference influences the 
photoresponse of the parallel-coupled QDs. 

The rest of the paper is organized as follows. First, 
in the following section, we establish the quantum rate 
equations for this system in the presence of a mag- 
netic field by employing the nonequilibrium Green's 
functioniiii^ Then in the section III, we calculate the 
current as a function of magnetic fluxes, and study the 
quantum dynamics of this system. The spectrum inves- 
tigation in RefEl pointed out that increasing strength of 
the additional bridges causes the total localization of the 
antibonding state due to the perfect destructive interfer- 
ence, and as a consequence the transport characteristic of 
the device reduces approximately to a single QD. In this 
paper we restrict our interest to the regime where there 
are two distinctly resolved peaks in the density of states 
spectra and both of the bonding and antibonding states 
can contribute to transport. Our numerical results show 
that the current in this regime still keeps the oscillation 
behavior with magnetic flux but the period changes from 
27r to 477 due to the interdot coupling. The temporal in- 
vestigation of the electron occupation probabilities in the 
two QDs shows that the conventional oscillation behavior 
in a two-level system can be destroyed by the additional 
bridges connecting the two QDs and two leads, and it 
can be recovered by applying a nonzero magnetic flux. 
Also in this section we study in detail the photoresponse 
of the system subject to a small irradiation and predict 
novel enclosed magnetic flux-controlled photon-assisted 
peaks in tunneling current, which can be attributed to 
the interference between the two pathways electrons go- 
ing through in this system. Finally, a summary is given 
in Section IV. 



II. MODEL AND FORMULATION 

We consider the parallel-coupled interacting QDs inter- 
ferometer connected to two normal leads as depicted in 
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Fig. 1. Only one bare energy level in each dot is involved 
in transport. The intradot electron-electron Coulomb in- 
teractions are assumed to be infinite but the interdot 
interaction U is finite. Namely, the state of two electrons 
occupied in the same QD is forbidden but two electrons 
dwelling in different QDs is permitted. 
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FIG. 1: Parallel-coupled quantum dots Aharonov-Bohm in- 
terferometer. 

For the sake of simplicity, we abandon the spin be- 
cause transport through this system is spin independent. 
Therefore, the available states and the corresponding en- 
ergies for the interferometer with two embedded QDs are: 
(1) the whole system is empty, |0)i|0)2, and the energy 
is zero; (2) the first QD is singly occupied, |l)i|0)2, and 
the energy is ei; (3) the second QD is singly occupied, 
|0)i|l)2, and the energy is €2; (4) both of the QDs are 
singly occupied, |l)i|l)2, and the energy is ci + €2 + U. 
We assign these Dirac brackets as operators: the slave- 
boson operators = |0)i|0)2, rf^ = |l)i|l)2 and the 
pseudo-fermion operators /| = |l)i|0)2, fl = |0)i|l)2. 
And obviously the explicit (anti)communicators of these 
auxiliary particles are»i^ 



poo = 



Pdd = 

Pl2 = 



where the statistical expectations involve the Fourier 
transformations of the time-diagonal parts of the ma- 
trix correlation functions in 2 x 2 space [G'^ j^{t,t')\ij = 



eS = e/j = et ^ ^ = de^ = df] = 0, (1) 

in association with the completeness relation 

e^e^dU^ flf^+ flf^=l. (2) 

The density matrix elements are expressed as poo = 
|0)i|0)22(0|i(0| = ete, pn - |l)i|0)22(0|i(l| = /|A, 

P22 = |0)l|l)22(l|l(0| - /2V2, Pdd = |l>l|l)22(l|l(l| = 

d'l'd, and pi2 — |0)i|l)22(0|i(l| = flfi- In terms of these 
slave particles operators, the Hamiltonian for this system 
can be written as 

rj.k 

+Wlf^+flf,) + {2ea + U)dU 

+ 5][14ie'^/4c^,(et/i+/td)+H.c.] 

k 

+ E[^^2e-'^/^ci,(et/2 + fid) + H.c] 
+ £[l^me-'^/'4,(et/^ + /trf) +H.C.] 

k 

+ E[^«2e'^/^4,(et/2 + fid) + H.c], (3) 
fc 

where c]^^, (c]^^,) are the creation (annihilation) operators 
for electrons with moment fc, and energy e^^ in the lead 
rj {r] = L, R). Vr/j denotes the hopping matrix element 
between the dot and the lead and ip = 27r$/$o accounts 
for the enclosed magnetic flux inside the AB interferom- 
eter ($0 — h/e is the flux quantum), tc is the interdot 
hopping coupling. 

We evaluate the statistical expectations of the rate of 
change of the density matrix elements pij with the Hamil- 
tonian and modified quantization Eq. (QJ. After te- 
dious but straightforward calculations, we obtain 



(4) 



(5) 
(6) 



(7) 



^(4(t')etW/,W), [Gtkit,t%, ^ ^{cl{t')f!mt)), 
[G<Jt,t%, EE z(/j(t')e(t')c,fe(<)), and [G<Jt,t%, ee 



m,e^e]) = i-^ {[VtG<,(t,0]ii - [G<,(t,i)Ve]n + [VtG<,(t, t)]22 - [G<,(t, i)Ve]22} , 

k 

m,flf,]) = i-^ {[G<,(t,<)V,],, - [VlG<Jt,t)]u + [vlG<,{t,t)U - [G<,it,t)Y,U} 

k 

+itc{Pu - Pli), 

m,dU]) = {[G<,(t,<)Vrf]n - [VlG<Jt,t)]n + [G<,(t, i)V,]22 - [V\G<Jt,t)]22}, 

k 

m,flf,]) = i-^ {[G<,(t,<)V,]i2 - [VlG<Jt,t)],2 + [VlG<,{t,t)U [G<,{t,t)V,U} 



+i{e2 - ei)pi2 + itcipu - P22), 



r 
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i{d^t')fj{t')c^f.{t)). With the help of the Langreth an- 
alytic continuation rulesji^ we can relate these hybrid 
Green's functions to the dressed Green's functions of the 
central region: 



+G<^{t,h)Vl^^gl{h,t% (8) 



where Ve and are two 2x2 niatrixs of the hopping 
elements defined by 



(9) 



and 



are the exact 



Green's functions of the ith lead without coupling to the 
device. These retarded (advanced) and lesser (greater) 
GFs for the central region are defined as: G^Ji^\t,t') = 



and G>^ 

and Oj = fjd if o 



±^9{±tTmO,it),0]it')}), G<^it,t')^^{0]{t')0,it)) 
{t, t') = -i{0.^ {t)0]{t')) with Oj ^ e^fj iio = e 
d. 

In the following derivation we perform a "gradient ex- 
pansion" of Eq. ((SJ, which is first introduced by Davies et 
al. to get the rate equation for resonant tunneling in the 
sequential regime»i& First define center-of-mass and rela- 
tive times hyT={t + t')/2 and t = {t - t')/2. Then we 
assume that functions vary rapidly in the relative time t 
but only slowly in the center-of-mass time T. Finally, we 
take a Fourier transform from t to uj, and the GFs G(t, t') 
in Eq. (jHJ becomes G{uj^T). According to Ref. |l^ the 
lowest-order gradient expansion is a good approximation 
for sequential resonant tunneling. Therefore, we just re- 
main the first term in the gradient expansion of the GFs 
G(w, T) in Eq. © and substitute these GFs and the hop- 
ping matrixs Ve/d into Eqs. Q-IZll. It is noted that the 
equal time in Eqs. (©-0 means t = or an integral over 
all UJ. Under the weak coupling (dot-lead tunneling and 
interdot hopping) assumption and slowly varying in time 
T, the GFs in the isolated two QDs system can be ex- 
pressed in terms of spectrum representation j Insert- 
ing these GFs into the Fourier forms of Eqs. (jHl), we then 
obtain the final quantum rate equations in the sequential 
tunneling regime. Because our interesting is focused on 
studying quantum dynamics and photo response of this 
interferometer at zero temperature and large bias volt- 
age, we will not intend to give the general expressions 
but the interesting readers could refer to our recent pa- 
per Ref. [T3 . Finally, at zero temperature and large bias 
voltage, the quantum rate equations are written as: 



+ [rm2e-''^/Vi2+H.c.], (10) 

Pll = ^LlP00+'^R2Pdd- (rRl+'^L2)Pll+itc{pi2- P2l) 

-[^(rm2e-^^/' + fL,2e'''/')pi2 + H.c], (11) 

P22 = rL2P00 + ^RlPdd - (rfl2 + rLl)p22 + itc{P21 ~ P12) 

-[^(rm2e~^^/' + rLi2e'^/2)pi2 + H.c], (12) 

Pdd = Ti^pii + riip22 - (r^i -I- TR2)pdd 

+ [rLi2e'^/V2 + H.c.], (13) 
Pi2 = «(e2 - ei)pi2 + »tc(pii - P22) + rLi2e""^/^poo 

+rR12e"^/^Pdd - ^(Ffll + Tr^ + f ii + fi2)pi2 

-^(rm2e^'^/' + rii2e-^'^/2)(pn + P22), (14) 



P22 
r . 



Poo 



FfllPll -I- Tr2P22 — (Fil + ri2)P00 



along with the normalization relation poo + pii ^ 
Pdd = 1 and P21 = p*2, with the definitions 
2^Sfc |K;iP<5(ti' — e-nk) denoting the strength of couphng 
between the ith QD level and the lead ry. Namely, Tn 
{TRi) here describes the tunneling rate of electrons in 
to (out from) the ith QD when the other QD is empty. 
On the contrary, Tn (F^) describes the tunneling rate 
of electrons in to (out from) the zth QDs, when the 
other QD is already occupied by an electron, revealing 
the modification of the corresponding rates due to the 
Coulomb repulsion between the two QDs. The interfer- 
ence in tunneling events through the different pathways 
are explicitly described by and F^y for the singly- 
occupied channel and doubly-occupied channel, respec- 
tively, with the definitions F,,ij = 27r VrjiVrij5{uj~ er^k) ■ 
These tunneling parameters are taken as constant under 
the wide band limit. In addition, the contribution of 
two leads is indeed negative to the nondiagonal density 
matrix element's dynamic equation, leading to damp- 
ing of the quantum superposition. It is obvious that 
these damping terms are different from the seirs-coupled 

QDsiaiiiia 

Actually, the similar equations have already been de- 
veloped for this system by using other schemes^iiS 
Jiang and co-workersii applied the Gurvitz's wavefunc- 
tion methodpiS to derive the modified rate equations and 
studied the temporary dynamics. However, it should be 
noted that their equations are different from ours for the 
nondiagonal density matrix element Eq. (|14l) . Marquardt 
and Brudesi^ started from the von Neumann equation of 
the reduced density matrix and obtained the rate equa- 
tions at a finite temperature. They studied the dephas- 
ing in sequential tunneling due to electron-phonon inter- 
action for the similar device without interdot hopping. 
Their equations at zero temperature coincide with ours 
in absence of the interdot hopping. 

The particle current 1^ flowing from the lead rj to the 
interferometer can be evaluated from the rate of change 
of the electron number operator iV,,(t) = Yl,k '^lik('^)^rik('^) 
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of the lead rj 



,20 



(15) 

Ultimately, the current I^/ji can be expressed in terms 
of the GFs: 

Il/h = g^{GLHVt-V,G<,H 



/22' 



(16) 



Under the weak coupling approximation, it becomes 



rLi2{e"^^^Pi2+e- 



^icp/2 



P2l) 



(17) 



The final term in the above expression comes from the 
contribution of the interference between the upper 
lower pathways. 
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FIG. 2: The calculated stationary current 7 as a function of 
the magnetic flux (p for different interdot coupling t — 0, 0.01, 
0.1, and 1 with T = 0.2 and T' = O.IF. 



of TT. But the positions of the current valleys move from 
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III. CALCULATIONS AND DISCUSSION! 

A. Dc current and Quantum dynamics 

In this section, we first calculate the current thr( 
the parallel-coupled QDs in the presence of magnetic 
and then study the quantum dynamical behavior of 
system. The density of states and linear conductam 
this system without the Coulomb interaction have 
explored in detail in absence of the magnetic flu 
literatureAS, It is found that when the system cha 
from a configuration in series to a completely sym 
rical parallel geometry, the tunneling through the ; 
bonding state is totally suppressed due to the pe 
destructive quantum interference between the diffe 
pathways through the systemiS With this point of \ 
we limit our investigation on the case of an asymrr 
parallel configuration to guarantee that both the b 
ing and antibonding states have contribution to ti 
port, in order to demonstrate the effect of interfer 
on the photoresponse clearly. Here we assume the 
nehng rates Tli = Tr2 = Tli = Tbs = T = 0.2 
rL2 = Tni = rL2 = Tri = r < O.Sr. And we 

rti2 = ^Ri2 — Vrr'. 

Figure 2 displays the stationary current calcul 
from Eq. H17|l as a function of the renormalized mag] 
flux if for different interdot couplings under the cond: 
r' = O.ir. It is clear that in the absence of the inte 
coupling tc — 0, the current / exhibits periodic os( 
tion with a period of 2tt. The current peaks appet 
the phases of (2n + 1)tt {n is an integer number), 
the current nearly vanishes at the phases of 2mT. Tl 
the main feature of the conventional AB effect. Howi 
when the interdot coupling turns on, the periodicil 
the AB oscillation becomes iw. In the new AB oscilh 
pattern, the first current peak also locates at the p 





5 10 15 20 5 10 15 20 

time time 

FIG. 3: The calculated time evolutions of the electron- 
occupation probabilities in the two QDs pn (solid lines) 
and P22 (dashed lines), (a)-(d): Different interdot couplings 
tc = 0, 0.01, 0.1, and 1 with F'/F = 0.1; (e)-(g): Different 
ratios of the two tunneling rates F'/F — 0, 0.5, and 1 with 
tc = 1; (h): F'/F = 1 and (f ^ n. The time unit is 1/5F. 

We then study the quantum dynamical behavoir of the 
system. In calculation, we assume that the system ini- 
tially occupies the empty state poo = 1- Figure 3 shows 



5 



the time evolutions of the electron probability densities 
Pii and P22 for the states |l)i|0)2 and |0)i|l)2. In the 
absence of the interdot hopping tc — 0, the probabil- 
ity Pu is nearly equal to 1, while P22 is nearly af- 
ter a long time, because the higher "injection" rate to 
the first QD makes electrons occupy this QD, and the 
higher "escape" rate of the second QD leads to no elec- 
trons staying at all [Fig. 3(a)]. When we turn on the 
interdot hopping tc, the probability in the second QD 
P22 is largely enhanced even for very small interdot hop- 
ping tc = 0.01 as depicted in Fig. 3(b). It is also shown 
that these probabilities display temporal oscillations at 
small time, meaning that the electron vibrates back and 
forth between the two QD. And the oscillation period is 
shortened by raising the interdot hopping tc- This is the 
well-known small-t oscillation behavior in the two-level 
system: pn ~ cos^ (tct) and P22 ^ sin^ {tct). However, 
we will find out some new characteristics in this small-i 
oscillation for the parallel-coupled QDs. As indicated in 
Fig. 3(d)-(f), the raising tunnehng rate of the additional 
pathway F' gradually destroys the oscillations of the two 
probabilities due to the increasing quantum interference 
between the additional pathway and the original one. For 
the equal tunneling rates of the two pathways, pu and P22 
are of course equal and show no oscillations [Fig. 3(g)]. 
This is a consequence of perfect destructive quantum in- 
terference between the different pathways through the 
parallel-coupled QDs. Furthermore, it is already known 
that applying the magnetic field will change the scatter- 
ing phase shift in every QDs, and then vary the interfer- 
ence patterns. As a result, the totally vanishing small-t 
oscillations can be recovered to some extent by threading 
a magnetic flux, which can be observed in Fig. 3(h) for a 
normalized magnetic fiux ip = t:. 



B. Photoresponse 

In this subsection, we study the time-dependent tun- 
neling through the AB interferometer in the nonlin- 
ear regimes. We assume that a time-dependent oscil- 
lation signal is applied to the two interacting QDs, so 
that the bare energy detuning becomes time-dependent 
£2—61 = eo + ^ cos nt, where S is the amplitude and 
51 the frequency of the externally applied signal. In the 
following, we use the approach developed by Stoof and 
Nazaro\»i£ to investigate the limiting case of a slight am- 
plitude, S <^ n,T]^/jf , the linear photoresponse. 

For simplicity, we assume the interdot Coulomb inter- 
action U is infinite, whereas only one electron can be 
found inside the system, so pdd = and all F arc equal 
to 0. The quantum rate equations l|ll|) - (|14ll simplify to 



Pu = Tpoo - r'pii + **c(/0i2 - P21) 
~^[rm2e-^'^/Vi2 + H.c.], 

P22 = F'poO - rp22 + «<c(p21 - P12) 

-J[rm2e-^'^/V2 + H.c.], 



(18) 
(19) 



P12 = «(e2 - ei)pi2 + itc(/Oii - /322) + rLi2e *'^/^poo 

-ir;^i2e^^/'(Pii + P22) - ^(r' + r)pi2, 

(20) 

with P00 + P11+P22 = 1- Correspondingly, the stationary 
current simplifies as: 



/L/e = -(F + F')(l-pii -P22). 



(21) 



We rewrite the simplified quantum rate equations in ma- 
trix notation: 



dp 

at 



(r 



Co + ^ COS nt)p + c. 



(22) 



where p = (pii, P22, P12, P2i)^, c = 
[r,r',TLi2e-"^/^,rLi2e"^/T, and T, T, cq, and S 
are the matrix forms of tunneling rates, hopping 
between dots, time-independent and time-dependent 
energy differences corresponding to Eqs. (|18|I - H20|) . The 
stationary solution of these equations without irradiation 
is expressed as: 



JO) 



-(r + T + eo)"'c. 



(23) 



Under the condition of small oscillation amplitude, the 
time-dependent density matrix elements can be expanded 
as: 



.(2) 



(24) 



where p^^^^ and p*^^) are the positive (negative) fre- 
quency part of the first order correction and frequency- 
independent second order correction to the stationary so- 
lution p^^\ respectively. They are proportional to the 
small amplitude of the oscillating signal S. Substituting 
the above equation into the time-dependent rate equa- 
tions and expanding according to the perturbation pa- 
rameter S, therefore, we obtain 



1(2) 



1 



T + T + €QTi^iy Sp 



(25) 



l{T + T + eo)-'Sip^'+'>+p<-'-^), (26) 



I being the unit matrix. The first order correction p'-'^^' 
provides oscillatory terms and has no contribution to the 
dc current (time average current) Eq. H21|l . The remain- 
ing lowest order contribution of the oscillating signal to 
the dc current comes from the second order correction 
p(2)_ jt is called photocurrent Ip^: 



(2) 



V-(r + r')[pri 



(2)1 
- P22 ] 



(27) 



In Fig. 4 we plot the calculated photoresponses sta- 
tionary current as a function of eo and the irradia- 
tion frequency f2 for a given interdot hopping tc — 1- 
Fig. 4(a) displays the result without the additional path- 
way F' = 0, which is just the coupled QDs in series. So 
the characteristic of this figure is the same as Fig. 2 in 



FIG. 4: The normalized photocurrent of the paraUel-coupled QDs, as a function of the bare level difference eo/tc between the 
two QDs and the frequency Q/tc of the irradiation, (a): The rate F' = and the renormalized magnetic flux y = 0; (b)-(c): 
F'/F = 0.1, the magnetic fluxes ip = 0, n, and 2tt. 



ReflTol: 1) Two branches of resonant satellite peaks for 
eo and satisfying il^ — eQ + (positive branch for 
eo > 0, negative branch for eo < 0), i.e., resonant PAT 
occurs when the emission or absorption of one photon 
can match the renormalized energy difference i/e^ + 
of the two levels; 2) No satellite peaks appearance for 
the case of frequencies below 2tc. This is because the 
condition of resonant PAT is never satisfied due to the 
lower energy hfl of the applied irradiation. They are 
the two general conditions for occurrence of the reso- 
nant PAT. However, when we switch on the additional 



pathway with a small rate F'/F = 0.1, we can observe 
that the photocurrent behavior is significantly altered as 
depicted in Fig. 4(b)-(c), although the two main char- 
acteristics above mentioned keep unchanged. It is clear 
in Fig. 4(b) that the positive branch of the PAT peaks is 
enhanced but the negative branch is suppressed nearly to 
zero amplitude. We attribute this phenomenon to the in- 
terference between the two pathways electrons can travel 
through in this new configuration. The scattering phase 
shifts are different when electrons pass through the two 
QDs with different bare energy levels. Consequently, the 



positive energy spacing eo > leads to constructive quan- 
tum interference. On the contrary, the negative energy 
difference eo < induces destructive quantum interfer- 
ence. This is the reason of the new resonant PAT pattern 
shown in Fig. 4(b). This explanation can be further sub- 
stantiated by the fact that the application of the mag- 
netic fluxes will change the interference fashion, and thus 
modify the photocurrent behavior. In Fig. 4(c) the cal- 
culated photocurrent is plotted for a given magnetic flux 
if = TT. An amazing finding is that the photocurrent 
becomes very similar to the result of the series-coupled 
QDs, except with a reduced magnitude. Recalling that 
the periodicity of the AB oscillation is An in the pres- 
ence of nonzero interdot hopping as pointed out in the 
above subsection, it is easy to imagine that we will obtain 
the opposite photoresponse to the situation depicted in 
Fig. 4(b) if the renormalized magnetic flux is set to be 
27r. The numerical results are plotted in Fig. 4(d). It 
is obvious that the negative branch of the resonant PAT 
peaks rises but the positive one declines. 

In the following, we study in detail what actually hap- 
pens for the reduced satellite peaks in the photoresponse 
stationary current by analyzing the evolution of the reso- 
nant PAT peaks as a function of the energy difference cq 
for different ratios of rates F'/F and renormalized mag- 
netic fluxes at a given frequency of the applied irra- 
diation (it must be bigger than or equal to 2tc)- We 
plot the calculated results in Fig. 5. Figure 5(a) is for 
the frequency of f2 = 3tc and several rate ratios in the 
absence of magnetic field. We can observe more clearly 
that the two PAT peaks are always located at the points 
e± = ± y^H.^ — 4i^ irrespective of whether the additional 
pathway is switched on or not. As the rate of the ad- 
ditional pathway is rising, the peak located at positive 
position e+ increasingly heightens and the e_ peak de- 
creases. But the shapes of peaks remain similar with 
those of F' = 0, in which case the two peaks are entirely 
identical and have a Lorentzian line shape near the reso- 
nance point e± as revealed in Refll3, until a stronger rate 
of the additional pathway F'/F = 0.25. For this maxi- 
mum rate considered here, the photocurrent displays as a 
negative dip with very small amplitude instead of a peak. 
It can be approximated by a Fano line shape. This behav- 
ior has a similar interpretation with those of the linear 
conductance and density of states when the configura- 
tion changes from series to parallel.^ The interfering of 
two tunneling channels opens prominent perspectives for 
the Fano effect, i.e., asymmetric line shape of current 
at resonant point. In addition, the smaller frequency of 
irradiation leads to the appearance of the dip in pho- 
tocurrent easier (at a smaller rate F'). Therefore, the 
Fano line shape can be observed more clearly at = 2tc 
as shown in Fig. 5(c). 

The effect of the magnetic field on the photoresponse 
can be clearly observed from Fig. 5(b) and (c). When 
the renormalized magnetic fiux ip is equal to tt, the Fano 
line shape of the resonant PAT peak absolutely changes 
back to the conventional Lorentzian shape if the reso- 
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FIG. 5: Evolution of the two satellite peaks for different ratios 
F'/F and different renormalized magnetic fluxes ip. (a): The 
plots are for F'/F = 0, 0.05, 0.1, and 0.25, respectively, at 
ip — and Q — 3tc; (b): The plots are for (p/n = 0, 0.5, 1.0, 
1.5, and 2.0 at F'/F = 0.1 and = 3tc; (c): The plots are the 
same as (b) but for Q = 2tc. For the sake of comparison, the 
result of r' = is denoted as the thin line in (c) . 



nant condition is satisfied. Finally, we conclude that for 
the moderate ratios of rates F'/F and frequencies of irra- 
diation, the photoresponse of the considered system ex- 
presses itself from the PAT Lorentzian peak to the PAT 
Fano peak by tuning the magnetic field threaded inside 
this interferometer. 



IV. CONCLUSION 

In summary, we have presented the quantum rate equa- 
tions in the sequential tunneling regime by means of the 
nonequilibrium Green's function for a mesoscopic AB 
ring with two tunneling-coupled QDs embedded in the 
two arms. Employing this set of quantum rate equations, 
we calculated the AB oscillation current, the temporal 
evolution of the electron-occupation probabilities in the 
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two QDs, and the dc photocurrent as the photoresponse 
in the presence of a weak irradiation, at zero temperature 
and large bias voltage between the source and the drain. 

Our numerical studies show that the permission of hop- 
ping between the two dots changes the AB oscillation 
period to 47r in comparison to the conventional 27r pe- 
riod oscillation observed in the typical AB effect. On 
the other hand, we find that the small-i oscillation of 
the electron-occupation probabilities is also established 
in the two QDs, with a damping amplitude controlled by 
the asymmetry of the parallel configuration. When the 
configuration is completely symmetrical the small-t oscil- 
lation is totally destroyed. Interestingly, this oscillation 
behavior is reobtained if we vary the enclosed magnetic 
flux due to the interference effect. 

Finally, we have evaluated in detail the dc transport 
through the parallel-coupled QDs subject to a small ex- 
ternal harmonic irradiation. It is found that, like in the 
series-coupled QDs, the photocurrent of this system ex- 
hibits extra resonant peaks when the frequency of the ex- 
ternal signal matches the energy difference between the 
discrete states. Moreover, one branch of the PAT peaks 
in photocurrent is enhanced, while another branch is sup- 
pressed, which is dependent on the enclosed magnetic 
flux. This behavior is a consequence of quantum inter- 
ference between the different pathways electrons can pass 
through, and does not exist for series-coupled QDs, where 



only one pathway exists. For some appropriate rate ra- 
tios r'/F and frequencies of signal fi, the PAT peaks of 
photcurrent are composed of Lorentzian and Fano line 
shapes at the resonant points, respectively, which can 
also be controlled by applying magnetic fields. 

In this paper, we study the time-dependent resonant 
tunneling through double QDs in a parallel configura- 
tion. This device has received wide attention in theoret- 
ical and experimental investigations at present, because 
its tunability makes it possible for further applications 
in quantum computation and quantum information. We 
hope our theoretical results about the novel PAT proper- 
ties can stimulate experimental studies on this particular 
problem. 
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